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SUMMARY 

The caoillarv pumped loop (CPL) system is one of the most desirable devices to 
dissipate heart energy in the radiation environment of the Space Station providi g 

temperature variations of the system are shown and more investigations are suggeste 
for further improvement. 


INTRODUCTION 

A caoillarv pumped loop (CPL) system consists of an evaporator pump, a 
condenser, a subcooler a vapor transport line, a liquid transport “™ anda f “ggf ^ 
rth a possible starter pump (Ku and others, 1966; 

1990) The condenser structured as a concentric core of a double tube heat excnang 
dissipates heat energy to an evaporator, which is theannulussecrt'onof the d 
exchanoer so that the fluid in the core condenses while the fluid in the annulus section 
pvanorates A CPL system equipped with a condenser having a two-phase/two-pha 
double tube heat dissipation structure is to be simulated numerically. A schematic block 
diagram is shown in Fig.1 and illustrates the control devices o i 

The system has two different types of evaporator pumpseacho^ 
units of evaporator pumps, double two-phase heat exchanger (^IDBTPHX) sy^m, 
which has six DBTPHX units as shown in Fig.2, a subcooler next to the DOTPHX sy 
with a non-condensible gas collector, long liquid and vapor lines, a ® u ^ y ^ of ^ , 
devices consisting of valves and fluid meters and a reservoir wth a starter pump. ^ 
order to have a manageable simulation system, certain assumptions are made. ( ) 
Mid plate evaporator pump replaces a hybrid evaporator so that two identical ^dplate 
evaporator pumps are a set of evaporator pumps; (2) fan-shape inlet and ex,t se^'ons 
of the DBTPHX (Fig.2) are assumed to be an assembly of six straight tubes in 
level of elevation; (3) the non-condensible gas(NCG) collector in the sudcool ® r d ° 
affect the performance of the system; (4) the flow control devices 'or '^ sy^em 's r^ 
included in the model in question; and (5) the storagetankb^^^ 
with constant properties. Such assumptions result in the block diagram as snown 
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Fig.3. 


Starting from the end of the subcooier, a few key points are identified as J1001, 

J1002 J1005 along with the Line 101, Plen 999, evaporator, Line 401, TEMP2C 

condenser, TEMP3C evaporator, TEMP2C subcooler and R114 chiller. A code known 
as SINDA'85/FLUINT (Cullimore and others, 1989) is employed for the simulation scheme 
and similar notations to those of SINDA'85/FLUINT are used in Fig.3. Details of these 
components consisting of an entire simulation model are illustrated in the following 
sections. 


EVAPORATOR PUMP SYSTEM 

Cross-sections of the evaporator pump in two directions are shown in Figs.4a and 
4b. Fig.4c illustrates how the liquid flow evaporates. The cold plate evaporator 
pump(EVP) with twelve(12) EVP's, twelve isolators are shown in Fig.4d, and a simplified 
model is indicated within the cross-hatched area in Fig.4d with one unit of the heat-pump 
and an isolator attached to it. Fig.4e explains the cross-section of the circular heat-pump 
and attached plates. Based on this figure of the combined cross-section, an approximate 
method for the extended surface theory is used. Namely, at the mid point of the plate, 
the temperature is the maximum and the temperature at the tube is a fixed temperature. 
Therefore, an approximate fin efficiency is used for an average temperature(or 
conductance) as far as the plate is concerned. The radial cross-section of the EVP has 
40 internal grooves in the outer shell. From a header, liquid enters the isolator, 
permeates the porous layer, reaching the core of the EVP. Once the liquid reaches the 
groove surface through permeation, heat transfer from the outside causes vaporization 
of the liquid. Evaporated fluid is pushed to the grooves(Fig.4c) and to the vapor header 
and enters the vapor transport line. Cullimore (1989) successfully demonstrated a 
numerical model by using a MACRO command, CAPPMP. Therefore, his method is 
employed for the EVP system. The capacity of the EVP system is assumed to be 400 
watts. 


DOUBLE TWO-PHASE HEAT-EXCHANGER SYSTEM 

In the DBTPHX, the condenser is a inner circular tube coupled with an evaporator 
which is the annulus section of the DBTPHX. The inner tube has axial grooves internally 
and externally and porous material layers occupy the space next to the grooves in the 
annular section and the core section, respectively, so that liquid, from capillary action can 
permeate the grooves and the porous layers. 

Six DBTPHXs connected in parallel function as the condenser and are designated 
as TEMP2C. The layer of porous material, Porex, enhances condensation in the internal 
grooves. Six DBTPHX evaporators connected in parallel are designated as TEMP3C. 
The layer of Porex directs the liquid flow in one desired direction. Porex has a 
permeability of 2.3x1 0E- 13. 


258 



From the inlet of TEMP2C a vapor enters i internal graves rol -the 

condenses on the surface of the porous layer and permeates ' * SbtpHX the 

the core cavitv and leaves the TEMP2C from the core section of the DBTPH)e From tne 

into the grooves, leaving TEMP3C from the annulus From ,hls exrt ’ P 088 

for a mixture of saturated liquid and vapor to leave TEMP3C. 

Figs 5a and 5b show an axial cross-section and a fadial ^oss-serton 
resoectivelv Six of the fan-shaped inlet and exit sections are represented by _a sing e 

consider that the quality of condenser fluid changes from unity to zero, while the qua y 
o° evaporator fluid changes from zero to unity For simplrcny the vanahons of fluid 
™.aiitw in TFMP2C and TEMP3C are assumed to be linear as shown in Fig.6. To t e 
^ a nnint nf TFMP2C and TEMP3C the qualities remain unchanged then change linear y 
to' th^location^f rom the^mid^poin? of the' DBTPHX. In the present stud^the len^ of 
the DBTPHX is divided into ten segments. Accordingly, the quality will be assigne 
TEMP2C and TEMP3C. 

SUBCOOLER AND R114 CHILLER 

ThP subcooler of the CPL has a complicated structure for a vapor trap to cope with 
non-conSensible gas ?n the subcooler. Considering that the amount of noncondensib e 
oas irrelatively small in comparison with the flow rate, the vapor trap is not modeled m 
theDresemstudy In the location of the trap, a flat plate is placed Similarly to the 
evaporator pump^late, the half of the plate with the attached subcooler and the RII A 
chiller is considered to be like a fin having the minimum temperature ^ the m'd point of 
the contacting area between the subcooler and the chiller (see Fig.7). The fin ' e^ae V 
!s aSumed to be 80 per cent. The temperature variation along the transverse diction 
is assumed to be rather small in comparison to that in the direction of the tube axis. For 
each leg of the subcooler and the chiller a uniform temperature is assumed to exist and 
the 180 degree bends connecting four legs are considered to be adiabatic. This s t 
handled with a MACRO command of SINDA’85/FLUINT. 

The subcooled liquid passes through the liquid transportation ^ 
considered adiabatic. A pressure drop through this passage is added to the system 
pressure loss. The liquid then enters the evaporator pump. 

The liquid reservoir is added to the system model as a plenum, 
properties as constant. A mixing process in the reservoir may not result in consta 
properties, nonetheless, it is assumed a steady state process. 
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INPUT FILE FOR THE COMPUTER PROGRAM 


An input file for SINDA’85/FLUINT was constructed based on the following 
Headers: 

1 . Header Options Data 

2. Header Node Data, Plate 

3. Header Conductor Data, Plate 

4. Header Control Data 

5. Header Source Data, Plate 

6. Header Flow Data,TEMP2C 

7. Header Flow Data,TEMP3C 

8. Header Flow Data.STLTS 

9. Header Flogic 0.TEMP2C 

10. Header Flogic 0.TEMP3C 

11. Header Flogic 1.TEMP2C 

12. Header Flogic 2.TEMP2C 

13. Header Variables 1, Plate 

14. Header Output Calls, TEMP2C 

15. Header Operation Data 
Build DBTPHX, Plate 

Build DBTPHX, TEMP2C,TEMP3C,STLTS 
Call Fastic 
Call Fwdbck 

16. Conditional Call for Restar 

17. Header Subroutine Data if any 

1 8. End of data 

Details of these headers are explained in the manual(Cullimoreand others, 1989). 
Diffusion submodels consist of the evaporator tube, its plate with the web, the internal 
tube, the external tube of the DBTPHX and its subcooler and chiller bodies with the 
plates. All of these submodels are represented by a diffusion model, PLATE. Node and 
Conductor Data sections have the initial temperature, capacitance and conductance of 
these submodels. The fluid submodels are TEMP2C, TEMP3C,and STLTL. TEMP2C 
represents the evaporator pump, condenser, its subcooler and the transportation lines. 
TEMP3C represents the buffer evaporator and STLTL represents the refrigerant chiller. 
T ransportation lines for liquid and vapor are adiabatic. The numbers for nodes, conductors, 
lumps and connections of the EVP s are in 200's and the remaining numbers used for the 
thermal and fluid submodels are listed in Table 1. 

Using the standard notation for SINDA'85/FLUINT, one unit of EVP is illustrated 
in Fig. 8 by using the notations of the code, whereby the macro command CAPPMP is 
applicable. 
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A symbolic diagram is drawn for the condenser/TEMP2C (he evaporator/TEMPSC, 
the subcooTerATEMPZC. and the H1 14 chiller in Fig. 9 with an aid of the standard 
♦ QiMnA'85/FLUINT for DBTPHX. The fluid transportation lines are modeled 

* a MACRO N co D mmand L UNl Th°e suboooler/Rt14 chiller is modeled by a MACRO 

command HX. 


Evaporator 

Vapor Transportation 

Condenser/TEMP2C 

Subcooler/TEMP2C 

Evaporator/TEMP3C 

R114 Chiller 

Storage Tank 


of Numbers 

Line 

100 to 199 
200 to 299 

Line 

400 to 499 
500 to 599 
600 to 699 
700 to 799 
800 to 899 
900 to 999 


RESULTS and DISCUSSIONS 

A start-up performance is run as the first trial with a prepared file shown above. 

With a surrounding temperature of 6.6 C, every component of 

same temperature except the storage tank which is at 29 C. A heat trad of 400 w 

applied to L cold plate EVP, and for a certain time penod, 

k simulated. Computation is divided into two time zones: the first 2 mmutes ana tne 
fniinwina 28 minutes In the 2 minute duration, the first step of the calculation wa 
Fastlc's^procedu re to obtain a stable initial condition for the 

a stable condition, the Fwdbck procedure takes over he oomixjajon. The Fa*csch« an 

lxhr.rpeCrets,ories a, those locations for to Mg " 

FWDBCK scheme. Some of the results by Neiswanger and others(1987) seem to snow 

that the trend of temperature rise is similar. 

At about 25 minutes, a steady state is reached: the highest temperature at Jhe 
evaporator plate is 39.7 C and the fluid temperature of the evaporator pump is 38.4 C. 
This tendency should be compared with that of an expenmental re y- 

Results of shorter segments of DBTPHX were not included and the accuracy of 
the current result has not been established with respect to the segment size.ome 
Darameters like the pressure, quality, heat transfer also have not been included^ Other 
operating conditions such as starting from a usual standard condition an ^ a combir^tw 
o? Fwdbck initially and STDSTL operations have not yet been tested with this input f e. 


261 



These analyses are ongoing at present and will continue until* satisfactory results are 
obtained within a reasonable error limit. 

In the process of computation, it was necessary to use the Fastic scheme initially 
with the linear relationship for the quality along the half length from the mid point in both 
TEMP2C and TEMP3C sides of DBTPHX. As a result, a stable result was obtained 
within a few number of iterations in the Fastic calculation. Thereby, the Fwdbck 
procedure was carried out further, eventually to lead to a steady state condition. 

CONCLUDING REMARKS 

A CPL system with a DBTPHX condenser was analyzed to simulate performance 
of the system with simplifying assumptions. The maximum temperature at the 
evaporator plate reached 39.7 C. Further parameter studies should be done for better 
results of temperature, pressure and quality distributions. No experimental results have 
been compared to the present results. That comparison will be critical for refining the 
input file. 

In general, improvements can be made by increasing the number of the segments 
for the DBTPHX for more accurate solution, because the process of Fwdbck gives stable 
solutions but does not provide accurate solutions. Similarly, the evaporator pump and the 
chiller legs should be divided into several segments and search for better solutions for 
performance simulation is desirable. Other modeling methods for the DBTPHX may be 
tested and their results should be compared with the result presented here. 

The effect of the neighboring isolators on the liquid flow, and precise heat 
convection coefficients at the evaporating and condensing surfaces in the grooves and 
porous layers are not available. Those values that come from the subroutine are the 
convection coefficients for the usual boiling and condensing conditions over a flat plate 
or cylinder(ASHRE,1 989; Chen, 1963). The heat transfer coefficients used in this paper 
may yield first order approximate performances, however, the average heat transfer 
coefficient obtained from the NASA experiment (Neiswanger and McIntosh, 1987; 
Cullimore, 1989) will be used for further testing. Finally, deleting the simplifying 
assumptions in making a working model will yield a satisfactory simulation. 
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Fig. 3 A S implied Block Diagram for the CPL System 
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Pig. 7 Sabcooler/TBMP2C,R114 Chiller end 
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Pig. 8 Block Diagram of CAPPMP for the Evaporator System 
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Fig. 10 Temperature-Histories at Evaporator, TEMP3C 
and Subcooler for the First Two Minutes 
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Fig. 11 Temperature-Histories at Evaporator, TEMP3C 
and Sucooler for 30 Minutes 
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